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Abstract 

Here we present a self-consistent quasi-particle model for quark-gluon plasma and apply it to explain 
the non-ideal behaviour seen in lattice simulations. The basic idea, borrowed from electrodynamic 
plasma, is that the gluons acquire mass as it propagates through plasma due to collective effects and is 
approximately equal to the plasma frequency. The statistical mechanics and thermodynamics of such 
a system is studied by treating it as an ideal gas of massive gluons. Since mass or plasma frequency 
depends on density, which itself is a thermodynamic quantity, the whole problem need to be solved 
self-consistently. 
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1 Introduction 



The quasiparticle model of quark-gluon plasma (qQGP) is used to study the thermal properties of quark- 
gluon plasma (QGP) ^ El El ^ El- QGP is a plasma of quarks and gluons which exhibits a collective 
behaviour because of the interactions, governed by quantum chromodynamics (QCD). At extreamly high 
temperature, the interaction may be very weak because of the asymptotic freedom and hence the plasma is 
almost ideal gas of quarks and gluons. However, near and above the critical Tc, the coupling constant is 
not weak and hence QGP may be non-ideal gas of quarks and gluons. Here Tc is the critical temperature of 
phase transition or cross over from hadrons to QGP. The lattice simulation of QCD at finite temperature [5] 
(LGT) and the elliptical flow observed in relativistic heavy ion collisions (RHICs) signal the existance of non- 
ideal QGP. The problem of confinement-deconfinement of QCD and the search for QGP in RHICs depend 
on the properties of QGP at temperature close to T^. Since the strong coupling constant oig is not very 
small near Tc, it may be difficult to solve QCD analytically. Many attempts in this line were not completely 
successful [TllSlin]. Therefore, several phenomenological models with two or more fitting parameters, like 
confinement models ^Uji strongly coupled plasma (SCQGP) 11 , strongly interacting plasma (sQGP) |12|. 
qQGP ^|21E101|S], liquid model P[3] and so on, were proposed. By adjusting parameters of the model 
one may fit the LGT results. Here we propose a new quasiparticle model which has only one adjustable 
parameter and is the revised version of the widely studied qQGP [SlITlini. 

2 Self-consistent qQGP: 

Let us recollect first various versions of qQGP models. It was first proposed by Peshier et. al. and with 
two adjustable parameters, they were were able to fit LGT results of |15| on gluon plasma. However, it fails 
to fit the more recent LGT results as shown in JT]. Furthermore, Gorenstein and Yang ^| pointed out that 
there is a thermodynamic (TD) inconsistency in this model and by imposing a stringent constraints, which 
they called TD consistency relation or the reformulation of statistical mechanics (SM), they tried to solve 
the problem. Recently, in a series of papers |0, we pointed out the reason for the TD inconsistency and 
revised their model using the standard SM in a TD consistent way. The revised model of qQGP, without TD 
consistency relation or the reformulation of SM, fits remarkable well many LGT results with a single system 
dependent parameter We also showed how the TD consistency relation comes out from our formulation. 
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The basic idea in our model is that we start from the energy density of quasiparticles instead of ideal gas 
pressure-partition function relation. Further, we assume that quasiparticlc acquire an additional temperature 
dependent mass and this thermal mass is equal to the plasma frequency, ujp, because of the collective effects 
of QGP. We used, as an approximation, the expression for ujp from perturbative QCD results and assumed 
the two-loop running coupling constant as function of temperature. So we just have one adjustable parameter 
related to QCD scale parameter which appears in as{T). In this letter, we further revise our qQGP model 
by considering a density dependent plasma frequency, instead of the perturbative QCD results which is valid 
at extreamly high temperature. It is motivated from a similar work in ultra-relativistic (e, e+ and 7) system 
|16[ll7j where one uses ultra-relativistic plasma frequency which depends on density and temperature. 

Here we consider a simplest and well studied system, gluon plasma (GV). It is a system of gluons 
interacting via QCD interaction at finite temperature, T. As we discussed earlier, it is difficult to solve 
it analytically since the coupling constant is not weak. Therefore, we consider a phenomcnological model 
where the thermal properties may be obtained by studying the thermal excitations of plasma modes. These 
thermal excitations are called quasi-gluons with the quantum numbers of gluons and with the thermal mass 
equal to the plasma frequency. Thus we have a gas of non-interacting or ideal quasi-gluons and following 
the standard SM , the logaritham of grand partition function or q-potential is given by 

00 

g = -^ln(l-e-'^^^) , (1) 

k=0 

where is the single particle energy of quasi-gluon, i.e, gluon with temperature dependent mass, given by, 

where k is the momentum and m is the thermal mass which is equal to the plasma frequency. f3 is defined 
as 1/T. Instead of using the QCD perturbative expression for plasma frequency as done in the earlier qQGP 
models 13 , we model it as 

cx - = ao a, - , (2) 

which is motivated from a similar work on electrodynamic plasma (EDP) I17| . The constant ao — ^ 
in electron-positron EDP plasma and here we fix it by demanding that ui^ ^ ^ as T — > 00, the QCD 
perturbative result, is related to as through the relation as = .g^/47r. Further, we use temperature 
dependent running coupling constant, as(r), which is motivated from lattice simulation of QCD n is 
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the density of gluons which is taken to be the same as the number of quasi-gluons in quasiparticle models. 
Again from standard SM, 



which may be rewritten as 



n=4^ dkP , , (3) 



/•OO -| 



where gf — 16, the degeneracy associated with the internal degrees of freedom. These equation need to be 
solved self-consistently because n which is to be determined is inside the integral through ujp. Note that we 
don't use the perturbative QCD expression for plasma frequency, which is appropriate at T oo, instead 
we calculate it self-consistently. Redefining the variables, the final equation to be solved self-consistently is, 

/•OC -1 -1 

f^^ dxx^ _^ ^c?flY.-K,{dlf,), (5) 

Jo ev^ - 1 1~[ ^ 



where 



and 



-27r2 



2 ^ Ztt^ n 



where K2 is the modified Bessel function. Above equation, Eq. Q, may be solved numerically to get fg 
and then other TD quantities like energy density, pressure etc. may be calculated. 
The energy density is given by, 

27r^ 7o ev + ■'^s - 1 

or 

00 

£ = /^^'E [(«/ffO'i^i(^/90 + 3(a/,0'^2(a/,/)]) , (7) 
^ 1=1 

in terms of modified Bessel functions K. Next, the pressure may be obtained from the TD relation, 

dP 

e^T—~P, (8) 

on integration, which is the procedure we follow here. In this analysis, we have neglected vacuum energy or 
zero point energy with the assumption that the whole thermal properties of giuon plasma may described by 
the quasi-gluon excitations. A new TD consistent qQGP model with the vacuum energy is presented in Ref. 
0. 
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Figure 1: Plots of P/T"^ (lower set of graphs) and e/T"^ (upper set of graphs) as a function of T/Tc from 
our model and lattice results (symbols) for gluon plasma. 

3 Results: 

In Fig. 1, we plotted energy density and pressure of GP from our model along with the LGT results |^. 
First we solve, self-consistently, the integral equation for the density, Eq. Q, for a given temperature and 
obtain fg{T). We have used 2-loop order running coupling constant, cks(T), which is similar to one used in 
LGT calculations, and is given by. 



6tt 



1 - 



3(153 - 19nf) ln(21n(r/AT)) 



(9) 



(33-2n/)ln(r/AT) V (33 - 2n/)2 ln(T/AT) 
where At is a parameter related to QCD scale parameter and Uf is the number of flavors which is zero 
in our case. Once we know fg(T), energy density and pressure are calculated using Eq. |(7J) and Eq. (jS)) 
respectively. We adjust only one parameter of our model, to = At/Tc, so that we get the best fit to the 
energy density of LGT results 0. We found that to — 0.83. We need one more integration constant, which 
is not a parameter, to evaluate pressure which we fix to LGT data at T = Tc. We may as well fix it at 
T — oo such that P/Psb — 1, where Psb is the ideal gas pressure, but it is difficult to solve numerically. 

In Fig. 2, we plotted fg{T), as{T) and u!p/{g'^ T^) as a function of T/Tc. The running coupling constant 
as{T) increases rapidly as T ^ Tc and oj^/ {g^ T^) is small near T — Tc and asymptotically approaches 1/3, 
the QCD perturbative result. 
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Figure 2: Plots of (T), ^^(T) and T^) as a function of T/T, 



4 Conclusions: 



We presented a new quasiparticle model for gluon plasma where gluons acquire mass, approximately equal 
to the plasma frequency, due to collective effects. Since ujp depends on density, which is one of the ther- 
modynamic quantity to be determined from statistical mechanics, we solved the problem self-consistently. 
Using this result, energy density and pressure were evaluated and by adjusting a single parameter, related 
to QCD scale parameter, we got a remarkable good fits to LGT results. Further extension of the model to 
QGP with quarks may be interesting, but not as simple as GP. 
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